Mitogenomes do not substantially improve phylogenetic resolution in a young non-model adaptive radiation of freshwater gastropods.

BACKGROUND: Species flocks in ancient lakes, and particularly those arising from adaptive radiation, make up the bulk of overall taxonomic and morphological diversity in these insular ecosystems. For these mostly young species assemblages, classical mitochondrial barcoding markers have so far been key to disentangle interspecific relationships. However, with the rise and further development of next-generation sequencing (NGS) methods and mapping tools, genome-wide data have become an increasingly important source of information even for non-model groups. RESULTS: Here, we provide, for the first time, a comprehensive mitogenome dataset of freshwater gastropods endemic to Sulawesi and thus of an ancient lake invertebrate species flock in general. We applied low-coverage whole-genome sequencing for a total of 78 individuals including 27 out of the 28 Tylomelania morphospecies from the Malili lake system as well as selected representatives from Lake Poso and adjacent catchments. Our aim was to assess whether mitogenomes considerably contribute to the phylogenetic resolution within this young species flock. Interestingly, we identified a high number of variable and parsimony-informative sites across the other 'non-traditional' mitochondrial loci. However, although the overall support was very high, the topology obtained was largely congruent with previously published single-locus phylogenies. Several clades remained unresolved and a large number of species was recovered polyphyletic, indicative of both rapid diversification and mitochondrial introgression. CONCLUSIONS: This once again illustrates that, despite the higher number of characters available, mitogenomes behave like a single locus and thus can only make a limited contribution to resolving species boundaries, particularly when introgression events are involved.

A trimeric glycosylated GH45 cellulase from the red abalone (Haliotis rufescens) exhibits endo and exoactivity.

The red abalone (Haliotis rufescens) represents North America's most important aquaculture species. Its hepatopancreas is rich in cellulases and other polysaccharide-degrading enzymes, which provide it the remarkable ability to digest cellulose-rich macroalgae; nevertheless, its cellulolytic systems are poorly explored. This manuscript describes some functional and structural properties of an endogenous trimeric glycosylated endoglucanase from H. rufescens. The purified enzyme showed a molecular mass of 23.4 kDa determined by MALDI-TOF mass spectrometry, which behaved as a homotrimer in gel filtration chromatography and zymograms. According to the periodic acid-Schiff reagent staining, detecting sugar moieties in SDS-PAGE gel confirmed that abalone cellulase is a glycoprotein. Hydrolysis of cello-oligosaccharides and p-nitrophenyl-ß-D-glucopyranosides confirmed its endo/exoactivity. A maximum enzyme activity toward 0.5% (w/v) carboxymethylcellulose of 53.9 ± 1.0 U/mg was achieved at 45°C and pH 6.0. We elucidated the abalone cellulase primary structure using proteases and mass spectrometry methods. Based on these results and using a bioinformatic approach, we identified the gene encoding this enzyme and deduced its full-length amino acid sequence; the mature protein comprised 177 residues with a calculated molecular mass of 19.1 kDa and, according to sequence similarity, it was classified into the glycosyl-hydrolase family 45 subfamily B. An AlphaFold theoretical model and docking simulations with cellopentaose confirmed that abalone cellulase is a ß-sheet rich protein, as also observed by circular dichroism experiments, with conserved catalytic residues: Asp26, Asn109, and Asp134. Interestingly, the AlphaFold-Multimer analysis indicated a trimeric assembly for abalone cellulase, which supported our experimental findings. The discovery and characterization of these enzymes may contribute to developing efficient cellulose bioconversion processes for biofuels and sustainable bioproducts.

Genome-Wide Identification and Transcriptome Analysis of the Hsp70 Gene Family in Monodonta labio Reveals Its Role in Response to Nanoplastics Stress.

For marine invertebrates, the disruption of organismal physiology and behavior by nanoplastics (NPs) has been extensively reported. Heat shock proteins (Hsps) are important for redundant protein breakdown, environmental changes, and intracellular protein transport. An exhaustive identification of Hsp70 genes and an experiment where different concentrations of NPs were stressed were performed to study how Hsp70 genes respond to NPs stress in Monodonta labio. Our results identified 15 members of Hsp70 within the genome of M. labio and provided insights into their responses to different concentrations of acute NP stress. Phylogenetic analyses revealed extensive amplification of the Hsp70 genes from the Hsc70 subfamily, with gene duplication events. As a result of NP stress, five of fifteen genes showed significant upregulation or downregulation. Three Hsp70 genes were highly expressed at an NP concentration of 0.1 mg/L, and no genes were downregulated. At 10 mg/L, they showed significant upregulation of two genes and significant downregulation of two genes. At 1 mg/L treatment, three genes were significantly downregulated, and no genes were significantly upregulated. Moreover, a purifying selection was revealed using a selection test conducted on duplicate gene pairs, indicating functional redundancy. This work is the first thorough examination of the Hsp70s in Archaeogastropoda. The findings improve knowledge of Hsp70s in molluscan adaptation to NP stress and intertidal living and offer essential data for the biological study of M. labio.

Characterization of the complete mitochondrial genome of Desmaulus extinctorium (Littorinimorpha, Calyptraeoidea, Calyptraeidae) and molecular phylogeny of Littorinimorpha.

For the purpose of determining the placement of Calyptraeidae within the Littorinimorpha, we hereby furnish a thorough analysis of the mitochondrial genome (mitogenome) sequence of Desmaulus extinctorium. This mitogenome spans 16,605 base pairs and encompasses the entire set of 37 genes, including 13 PCGs, 22 tRNAs and two rRNAs, with an evident AT bias. Notably, tRNASer1 and tRNASer2 lack dihydrouracil (DHU) arms, resulting in an inability to form a secondary structure. Similarly, tRNAAla lacks a TΨC arm, rendering it incapable of forming a secondary structure. In contrast, the remaining tRNAs demonstrate a characteristic secondary structure reminiscent of a cloverleaf. A comparison with ancestral gastropods reveals distinct differences in three gene clusters (or genes), encompassing 15 tRNAs and eight PCGs. Notably, inversions and translocations represent the major types of rearrangements observed in D. extinctorium. Phylogenetic analysis demonstrates robust support for a monophyletic grouping of all Littorinimorpha species, with D. extinctorium representing a distinct Calyptraeoidea clade. In summary, this investigation provides the first complete mitochondrial dataset for a species of the Calyptraeidae, thus providing novel insights into the phylogenetic relationships within the Littorinimorpha.

The genome of the rayed Mediterranean limpet Patella caerulea (Linnaeus, 1758).

Patella caerulea (Linnaeus, 1758) is a mollusc limpet species of the class Gastropoda. Endemic to the Mediterranean Sea, it is considered a keystone species due to its primary role in structuring and regulating the ecological balance of tidal and subtidal habitats. It is currently being used as a bioindicator to assess the environmental quality of coastal marine waters and as a model species to understand adaptation to ocean acidification. Here, we provide a high-quality reference genome assembly and annotation for P. caerulea. We generated ∼30 Gb of Pacific Biosciences high-fidelity data from a single individual and provide a final 749.8 Mb assembly containing 62 contigs, including the mitochondrial genome (14,938 bp). With an N50 of 48.8 Mb and 98% of the assembly contained in the 18 largest contigs, this assembly is near chromosome-scale. Benchmarking Universal Single-Copy Orthologs scores were high (Mollusca, 87.8% complete; Metazoa, 97.2% complete) and similar to metrics observed for other chromosome-level Patella genomes, highlighting a possible bias in the Mollusca database for Patellids. We generated transcriptomic Illumina data from a second individual collected at the same locality and used it together with protein evidence to annotate the genome. A total of 23,938 protein-coding gene models were found. By comparing this annotation with other published Patella annotations, we found that the distribution and median values of exon and gene lengths was comparable with other Patella species despite different annotation approaches. The present high-quality P. caerulea reference genome, available on GenBank (BioProject: PRJNA1045377; assembly: GCA_036850965.1), is an important resource for future ecological and evolutionary studies.

Molecular and structural analysis of Hdh-MIRP3 and its impact on reproductive regulation in female Pacific abalone, Haliotis discus hannai.

Molluscan insulin-related peptides (MIRP) play a crucial role in various biological processes, including reproduction and larval development in mollusk species. To investigate the involvement of MIRP in the ovarian development of Pacific abalone (Haliotis discus hannai), the Hdh-MIRP3 was cloned from cerebral ganglion (CG). Hdh-MIRP3 cDNA was 993 bp long, encoded a 13.22 kDa peptide, comprising 118 amino acids. Fluorescence in situ hybridization confirmed the localization of Hdh-MIRP3 in the CG and ovary. Molecular docking revealed that Hdh-MIRP3 binds to the N-terminal region of Hdh-IRP-R. Tissue expression analysis showed the highest Hdh-MIRP3 expression in the CG, followed by ovarian tissue. Hdh-MIRP3 expression was significantly upregulated in the CG and ovary during the ripe stage of seasonal ovarian development and in effective accumulative temperature conditioned abalone. Furthermore, siRNA silencing of Hdh-MIRP3 significantly downregulated the expression of four reproduction-related genes, including Hdh-GnRH, Hdh-GnRH-R, Hdh-IRP-R, and Hdh-VTG in both the CG and ovary, and Hdh-MIRP3 as well. These results indicate that Hdh-MIRP3 acts as a regulator of ovarian development in Pacific abalone. Additionally, expression analysis indicated that Hdh-MIRP3 plays a role in embryonic and larval development. Overall, the present findings elucidate the role of Hdh-MIRP3 in reproductive development in female Pacific abalone.

Toxic effects of combined exposure to homoyessotoxin and nitrite on the survival, antioxidative responses, and apoptosis of the abalone Haliotis discus hannai.

Homoyessotoxin (homo-YTX) and nitrite (NO2-N), released during harmful dinoflagellate cell lysis adversely affect abalones. However, their toxicity mechanisms in shellfish remain unclear. This study investigated the economic abalone species Haliotis discus hannai exposed to varying concentrations of homo-YTX (0, 2, 5, and 10 µg L-1) and NO2-N (0, 3, and 6 mg L-1) on the basis of their 12 h LC50 values (5.05 µg L-1 and 4.25 mg L-1, respectively) and the environmentally relevant dissolved concentrations during severe dinoflagellate blooms, including mixtures. The test abalones were exposed to homo-YTX and NO2-N for 12 h. The mortality rate (D), reactive oxygen species (ROS) levels, antioxidant defense capabilities, and expression levels of antioxidant-related, Hsp-related, and apoptosis-related genes in abalone gills were assessed. Results showed that the combined exposure to homo-YTX and NO2-N increased the D and ROS levels and upregulated B-cell lymphoma-2 (BCL2)-associated X (BAX) and caspase3 (CASP3) expression levels while reducing glutathione peroxidase (GPx) activity and GPx, CuZnSOD, and BCL2 expression levels. High concentrations of homo-YTX (10 µg L-1) and NO2-N (6 mg L-1) solutions and the combinations of these toxicants inhibited the activities of superoxide dismutase (SOD) and catalase (CAT) and downregulated the expression levels of MnSOD, CAT, Hsp70, and Hsp90. The ROS levels were negatively correlated with the activities of SOD, CAT, and GPx and the expression levels of MnSOD, CuZnSOD, CAT, GPx, Hsp70, Hsp90, and BCL2. These results suggest that homo-YTX, in conjunction with NO2-N, induces oxidative stress, disrupts antioxidant defense systems, and triggers caspase-dependent apoptosis in the gills of abalone. ROS-mediated antioxidative and heat-shock responses and apoptosis emerge as potential toxicity mechanisms affecting the survival of H. discus hannai due to homo-YTX and NO2-N exposure.

Chromosome-level genome assembly of Oncomelania hupensis: the intermediate snail host of Schistosoma japonicum.

BACKGROUND: Schistosoma japonicum is a parasitic flatworm that causes human schistosomiasis, which is a significant cause of morbidity in China, the Philippines and Indonesia. Oncomelania hupensis (Gastropoda: Pomatiopsidae) is the unique intermediate host of S. japonicum. A complete genome sequence of O. hupensis will enable the fundamental understanding of snail biology as well as its co-evolution with the S. japonicum parasite. Assembling a high-quality reference genome of O. hupehensis will provide data for further research on the snail biology and controlling the spread of S. japonicum. METHODS: The draft genome was de novo assembly using the long-read sequencing technology (PacBio Sequel II) and corrected with Illumina sequencing data. Then, using Hi-C sequencing data, the genome was assembled at the chromosomal level. CAFE was used to do analysis of contraction and expansion of the gene family and CodeML module in PAML was used for positive selection analysis in protein coding sequences. RESULTS: A total length of 1.46 Gb high-quality O. hupensis genome with 17 unique full-length chromosomes (2n = 34) of the individual including a contig N50 of 1.35 Mb and a scaffold N50 of 75.08 Mb. Additionally, 95.03% of these contig sequences were anchored in 17 chromosomes. After scanning the assembled genome, a total of 30,604 protein-coding genes were predicted. Among them, 86.67% were functionally annotated. Further phylogenetic analysis revealed that O. hupensis was separated from a common ancestor of Pomacea canaliculata and Bellamya purificata approximately 170 million years ago. Comparing the genome of O. hupensis with its most recent common ancestor, it showed 266 significantly expanded and 58 significantly contracted gene families (P < 0.05). Functional enrichment of the expanded gene families indicated that they were mainly involved with intracellular, DNA-mediated transposition, DNA integration and transposase activity. CONCLUSIONS: Integrated use of multiple sequencing technologies, we have successfully constructed the genome at the chromosomal-level of O. hupensis. These data will not only provide the compressive genomic information, but also benefit future work on population genetics of this snail as well as evolutional studies between S. japonicum and the snail host.

Cloning and characterization of heat shock transcription factor 1 and its functional role for Hsp70 production in the sea slug Onchidium reevesii.

To investigate the regulatory role of heat shock transcription factor 1 of sea slug Onchidium reevesii (OrHSF1) on Hsp70 expression in the sea slug under stress , the OrHSF1 gene was cloned and bioinformatics analysis was performed, then the gene and protein expressions by RNA interference (RNAi) mediated knockdown of OrHSF1 expression were measured to clarify the regulatory relationship between OrHSF1 and Hsp70 under low-frequency noise (LFN) stress. Our study was the first to clone a 1572 bp sequence of the OrHSF1 gene, with the sequence coding for amino acids (CDS) being 729 bp, encoding 243 amino acids. O. reevesii shared a close evolutionary relationship with mollusks such as the Aplysia californica. OrHSF1 gene is widely expressed in different tissues of sea slugs, with the highest expression in the intestine and the lowest in the reproductive glands. Furthermore, we used RNA interference (RNAi) as a tool to silence the OrHSF1 gene in the central nervous system (CNS) and the results indicated that gene silencing was occurring systematically in the CNS and the suppression of OrHSF1 expression by RNAi-mediated gene silencing altered the expression of Hsp70; besides, the expression trends of OrHSF1 gene and Hsp70 were consistent in the 3 and 5-day RNAi experiment. Moreover, in sea slugs injected with siHSF1 and exposed to LFN, the mRNA expression and protein expression of Hsp70 in the CNS were significantly decreased compared to the low-frequency noise group (P < 0.05). This study demonstrated that OrHSF1 regulates Hsp70 expression in marine mollusks under low-frequency noise, and HSF1-Hsp70 axis plays a key role in stress response.

A chromosome-level genome for the nudibranch gastropod Berghia stephanieae helps parse clade-specific gene expression in novel and conserved phenotypes.

BACKGROUND: How novel phenotypes originate from conserved genes, processes, and tissues remains a major question in biology. Research that sets out to answer this question often focuses on the conserved genes and processes involved, an approach that explicitly excludes the impact of genetic elements that may be classified as clade-specific, even though many of these genes are known to be important for many novel, or clade-restricted, phenotypes. This is especially true for understudied phyla such as mollusks, where limited genomic and functional biology resources for members of this phylum have long hindered assessments of genetic homology and function. To address this gap, we constructed a chromosome-level genome for the gastropod Berghia stephanieae (Valdés, 2005) to investigate the expression of clade-specific genes across both novel and conserved tissue types in this species. RESULTS: The final assembled and filtered Berghia genome is comparable to other high-quality mollusk genomes in terms of size (1.05 Gb) and number of predicted genes (24,960 genes) and is highly contiguous. The proportion of upregulated, clade-specific genes varied across tissues, but with no clear trend between the proportion of clade-specific genes and the novelty of the tissue. However, more complex tissue like the brain had the highest total number of upregulated, clade-specific genes, though the ratio of upregulated clade-specific genes to the total number of upregulated genes was low. CONCLUSIONS: Our results, when combined with previous research on the impact of novel genes on phenotypic evolution, highlight the fact that the complexity of the novel tissue or behavior, the type of novelty, and the developmental timing of evolutionary modifications will all influence how novel and conserved genes interact to generate diversity.

Cloning, characterization, and spatio-temporal expression patterns of HdhSPARC and its responses to multiple stressors.

SPARC is an extracellular Ca2+-binding, secreted glycoprotein that plays a dynamic role in the growth and development of organisms. This study aimed to describe the isolation, characterization, and expression analysis of HdhSPARC in Pacific abalone (Haliotis discus hannai) to infer its potential functional role. The isolated HdhSPARC was 1633 bp long, encoding a polypeptide of 284 amino acid residues. Structurally, the SPARC protein in abalone is comprised of three biological domains. However, the structure of this protein varied between vertebrates and invertebrates, as suggested by their distinct clustering patterns in phylogenetic analysis. In early development, HdhSPARC was variably expressed, and higher expression was found in veliger larvae. Moreover, HdhSPARC was highly expressed in juvenile abalone with rapid growth compared to their slower-growing counterparts. Among the testicular development stages, the growth stage exhibited higher HdhSPARC expression. HdhSPARC was also upregulated during muscle remodeling and shell biomineralization, as well as in response to different stressors such as heat shock, LPS, and H2O2 exposure. However, this gene was downregulated in Cd-exposed abalone. The present study first comprehensively characterized the HdhSPARC gene, and its spatio-temporal expressions were analyzed along with its responses to various stressors.

Insights into the genome of the 'Loco' Concholepas concholepas (Gastropoda: Muricidae) from low-coverage short-read sequencing: genome size, ploidy, transposable elements, nuclear RNA gene operon, mitochondrial genome, and phylogenetic placement in the family Muricidae.

BACKGROUND: The Peruvian 'chanque' or Chilean 'loco' Concholepas concholepas is an economically, ecologically, and culturally important muricid gastropod heavily exploited by artisanal fisheries in the temperate southeastern Pacific Ocean. In this study, we have profited from a set of bioinformatics tools to recover important biological information of C. concholepas from low-coverage short-read NGS datasets. Specifically, we calculated the size of the nuclear genome, ploidy, and estimated transposable elements content using an in silico k-mer approach, we discovered, annotated, and quantified those transposable elements, we assembled and annotated the 45S rDNA RNA operon and mitochondrial genome, and we confirmed the phylogenetic position of C. concholepas within the muricid subfamily Rapaninae based on translated protein coding genes. RESULTS: Using a k-mer approach, the haploid genome size estimated for the predicted diploid genome of C. concholepas varied between 1.83 Gbp (with kmer = 24) and 2.32 Gbp (with kmer = 36). Between half and two thirds of the nuclear genome of C. concholepas was composed of transposable elements. The most common transposable elements were classified as Long Interspersed Nuclear Elements and Short Interspersed Nuclear Elements, which were more abundant than DNA transposons, simple repeats, and Long Terminal Repeats. Less abundant repeat elements included Helitron mobile elements, 45S rRNA DNA, and Satellite DNA, among a few others.The 45S rRNA DNA operon of C. concholepas that encodes for the ssrRNA, 5.8S rRNA, and lsrRNA genes was assembled into a single contig 8,090 bp long. The assembled mitochondrial genome of C. concholepas is 15,449 bp long and encodes 13 protein coding genes, two ribosomal genes, and 22 transfer RNAs. CONCLUSION: The information gained by this study will inform the assembly of a high quality nuclear genome for C. concholepas and will support bioprospecting and biomonitoring using environmental DNA to advance development of conservation and management plans in this overexploited marine snail.

Phylogenomic analyses reveal a single deep-water colonisation in Patellogastropoda.

Patellogastropoda, the true limpets, is a major group of gastropods widely distributed in marine habitats from the intertidal to deep sea. Though important for understanding their evolutionary radiation, the phylogenetic relationships among the patellogastropod families have always been challenging to reconstruct, with contradictory results likely due to insufficient sampling. Here, we obtained mitogenomic and phylogenomic data (transcriptomic or genomic) from six species representing the three predominantly deep-water patellogastropod families: Lepetidae, Neolepetopsidae, and Pectinodontidae. By using various phylogenetic methods, we show that mitogenome phylogeny recovers monophyly of eight families in most of the trees, though the relationships among families remain contentious. Meanwhile, a more robust family-level topology consistent with morphology was achieved by phylogenomics. This also reveals that these mainly deep-water families are monophyletic, suggesting a single colonisation of the deep water around the Jurassic. We also found a lack of significant correlation between genome size and habitat depth, despite some deep-water species exhibiting larger genome sizes. Our phylogenomic tree provides a stable phylogenetic backbone for Patellogastropoda that includes seven of the nine recognized families and paves the way for future evolutionary analyses in this major group of molluscs.

Fifth mass extinction event triggered the diversification of the largest family of freshwater gastropods (Caenogastropoda: Truncatelloidea: Hydrobiidae).

The fifth mass extinction event (MEE) at the Cretaceous-Palaeogene (K-Pg) boundary 66 million years ago (Ma) led to massive species loss but also triggered the diversification of higher taxa. Five models have been proposed depending on whether this diversification occurred before, during or after the K-Pg boundary and the rate of species accumulation. While the effects of the K-Pg MEE on vertebrate evolution are relatively well understood, the impact on invertebrates, particularly in freshwater ecosystems, remains controversial. One example is the hyperdiverse Hydrobiidae-the most species-rich family of freshwater gastropods. Whereas some studies place its origin in the Jurassic or even Carboniferous, most fossil records postdate the K-Pg event. We therefore used robustly time-calibrated multi-locus phylogenies of >400 species representing >100 hydrobiid genera to unravel its evolutionary history and patterns of diversification. We found that the family started diversifying shortly after the K-Pg boundary (∼60 Ma; 95% highest posterior density 52-69 Ma). Lineage richness gradually increased to the present and phylogenetic diversity until ∼25 Ma. These findings suggest that diversification was not initially driven by ecological opportunity. Combining the two criteria of timing and rate of diversification, a soft-explosive diversification model of aquatic vertebrates best fits the patterns observed. We also show that most higher hydrobiid taxa (i.e. subfamilies) diversified from the Middle Oligocene to Middle Miocene (i.e. 12-28 Ma). Two of the 15 major clades delimited are described here as new subfamilies (i.e. Bullaregiinae n. subfam. and Pontobelgrandiellinae n. subfam.), whose members are restricted to subterranean waters. Our results are an important contribution to understanding how the fifth MEE has shaped evolution and patterns of biodiversity in continental aquatic systems. Given the high extinction risks faced by many hydrobiids today, they also emphasise the need to study the biodiversity of vulnerable ecosystems.

Exposure to bisphenol A and fiber-type microplastics induce oxidative stress and cell damage in disk abalone Haliotis discus hannai: Bioaccumulation and toxicity.

Along with environmental pollution caused by rapid economic development and industrialization, plastic waste is emerging as a global concern in relation to marine ecosystems and human health. Among the microplastics, fiber-type microfibers (MF) and bisphenol A (BPA), which are widely used as plasticizers, do not decompose well in the ocean, and tend to accumulate in organisms, generating an increased oxidative stress response. This study investigated the abalones' antioxidant and cell death responses following exposure to the environmental pollutants MF and BPA. Levels of malondialdehyde (MDA) and DNA damage increased over time, demonstrating the degree of lipid peroxidation and DNA damage in abalones exposed to individual and combined environmental conditions of MF and BPA. Compared to the single MF and BPA exposure groups, the combined exposure group showed a higher expression of antioxidant enzymes. A similar pattern was seen in the expression of the apoptosis enzyme caspase-3. Both MF and BPA caused oxidative stress and antioxidant enzymes were expressed to alleviate it, but it is believed that cell damage occurred because the stress level exceeded the allowed range.

Comparative genomic analysis of chemosensory-related gene families in gastropods.

Chemoreception is critical for the survival and reproduction of animals. Except for a reduced group of insects and chelicerates, the molecular identity of chemosensory proteins is poorly understood in invertebrates. Gastropoda is the extant mollusk class with the greatest species richness, including marine, freshwater, and terrestrial lineages, and likely, highly diverse chemoreception systems. Here, we performed a comprehensive comparative genome analysis taking advantage of the chromosome-level information of two Gastropoda species, one of which belongs to a lineage that underwent a whole genome duplication event. We identified thousands of previously uncharacterized chemosensory-related genes, the majority of them encoding G protein-coupled receptors (GPCR), mostly organized into clusters distributed across all chromosomes. We also detected gene families encoding degenerin epithelial sodium channels (DEG-ENaC), ionotropic receptors (IR), sensory neuron membrane proteins (SNMP), Niemann-Pick type C2 (NPC2) proteins, and lipocalins, although with a lower number of members. Our phylogenetic analysis of the GPCR gene family across protostomes revealed: (i) remarkable gene family expansions in Gastropoda; (ii) clades including members from all protostomes; and (iii) species-specific clades with a substantial number of receptors. For the first time, we provide new and valuable knowledge into the evolution of the chemosensory gene families in invertebrates other than arthropods.

Systematics and Distributions of Upper Bathyal Species in Bathyancistrolepis, a Deep-Sea Whelk Genus Endemic to the Northwest Pacific (Gastropoda: Buccinidae).

The deep-sea buccinid snail genus Bathyancistrolepis is redefined based on the reconstruction of a molecular phylogeny and morphological examination of shell and radular characters. This genus is distinguished from other genera of the subfamily Parancistrolepidinae with a combination of shell traits, including (1) a low spire, (2) sharp, carinate spiral cords or keels and (3) a long, curved siphonal canal, but not with a difference in radular morphology as suggested by previous authors. Three allopatric or parapatric species are recognized in the upper bathyal (447-2057 m) waters around Japan and Taiwan: B. tokoyodaensis from off Hokkaido to Sagami Bay in the Northwest Pacific, B. trochoidea off Kumano-nada to Miyazaki in the Northwest Pacific and along Nansei Islands in the East China Sea, and B. taiwanensis sp. nov. in the South China Sea. These species bear large paucispiral protoconchs that are indicative of benthic early development without a pelagic larval period, and hence low dispersal capability. Seafloor topography seems to have acted as a barrier for their dispersal; the range of B. tokoyodaensis supports the previous finding that Izu Peninsula delimits westward distribution of bathyal gastropod species of boreal origins.

Genomic DNA extraction optimization and validation for genome sequencing using the marine gastropod Kellet's whelk.

Next-generation sequencing technologies, such as Nanopore MinION, Illumina Hiseq and Novaseq, and PacBio Sequel II, hold immense potential for advancing genomic research on non-model organisms, including the vast majority of marine species. However, application of these technologies to marine invertebrate species is often impeded by challenges in extracting and purifying their genomic DNA due to high polysaccharide content and other secondary metabolites. In this study, we help resolve this issue by developing and testing DNA extraction protocols for Kellet's whelk (Kelletia kelletii), a subtidal gastropod with ecological and commercial importance, by comparing four DNA extraction methods commonly used in marine invertebrate studies. In our comparison of extraction methods, the Salting Out protocol was the least expensive, produced the highest DNA yields, produced consistent high DNA quality, and had low toxicity. We validated the protocol using an independent set of tissue samples, then applied it to extract high-molecular-weight (HMW) DNA from over three thousand Kellet's whelk tissue samples. The protocol demonstrated scalability and, with added clean-up, suitability for RAD-seq, GT-seq, as well as whole genome sequencing using both long read (ONT MinION) and short read (Illumina NovaSeq) sequencing platforms. Our findings offer a robust and versatile DNA extraction and clean-up protocol for supporting genomic research on non-model marine organisms, to help mediate the under-representation of invertebrates in genomic studies.

Viral Diversity in Samples of Freshwater Gastropods Benedictia baicalensis (Caenogastropoda: Benedictiidae) Revealed by Total RNA-Sequencing.

Previously, the main studies were focused on viruses that cause disease in commercial and farmed shellfish and cause damage to food enterprises (for example, Ostreavirusostreidmalaco1, Aurivirus haliotidmalaco1 and Aquabirnavirus tellinae). Advances in high-throughput sequencing technologies have extended the studies to natural populations of mollusks (and other invertebrates) as unexplored niches of viral diversity and possible sources of emerging diseases. These studies have revealed a huge diversity of mostly previously unknown viruses and filled gaps in the evolutionary history of viruses. In the present study, we estimated the viral diversity in samples of the Baikal endemic gastropod Benedictia baicalensis using metatranscriptomic analysis (total RNA-sequencing); we were able to identify a wide variety of RNA-containing viruses in four samples (pools) of mollusks collected at three stations of Lake Baikal. Most of the identified viral genomes (scaffolds) had only distant similarities to known viruses or (in most cases) to metagenome-assembled viral genomes from various natural samples (mollusks, crustaceans, insects and others) mainly from freshwater ecosystems. We were able to identify viruses similar to those previously identified in mollusks (in particular to the picornaviruses Biomphalaria virus 1 and Biomphalaria virus 3 from the freshwater gastropods); it is possible that picorna-like viruses (as well as a number of other identified viruses) are pathogenic for Baikal gastropods. Our results also suggested that Baikal mollusks, like other species, may bioaccumulate or serve as a reservoir for numerous viruses that infect a variety of organisms (including vertebrates).